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Velocity relaxation of a particle in a confined compressible fluid 
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(Dated: 7 March 2013) 

The velocity relaxation of an impulsively forced spherical particle in a fluid confined by two parallel plane 
walls is studied using a direct numerical simulation approach. During the relaxation process, the momentum 
of the particle is transmitted in the ambient fluid by viscous diffusion and sound wave propagation, and 
the fluid flow accompanied by each mechanism has a different character and affects the particle motion 
differently. Because of the bounding walls, viscous diffusion is hampered, and the accompanying shear flow 
is gradually diminished. However, the sound wave is repeatedly reflected and spreads diffusely. As a result, 
the particle motion is governed by the sound wave and backtracks differently in a bulk fluid. The time 
when the backtracking of the particle occurs changes non-monotonically with respect to the compressibility 
factor e — vjac and is minimized at the characteristic compressibility factor. This factor depends on the 
wall spacing, and the dependence is different at small and large wall spacing regions based on the different 
mechanisms causing the backtracking. 
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I. INTRODUCTION 

The dynamics of particles dispersed in a fluid flowing 
through a microchannel are important for many chemical 
engineering processes, such as membrane separation and 
microfluidics. The dynamics of fluids and dispersed parti- 
cles are significantly affected by bounding walls; thus, the 
transport properties of dispersions in some simple bound- 
ing geometries have been studied using analytical meth- 
ods^. While such studies were limited to a steady flow 
system, recently, the unsteady dynamics of a dispersed 
particle in a confined fluid have been investigatedS^i^. 

In a dispersion system, the momentum of a particle 
is propagated in the ambient fluid via two mechanisms: 
viscous diffusion and sound propagation. Each of these 
mechanisms is accompanied by fluid flow of different 
character, which affects the particle motion. Here, we 
consider the velocity relaxation of a particle after adding 
an impulsive force. In a bulk fluid, part of the particle 
momentum is transported by a sound wave to an infinite 
distance in time, and finally, the particle motion is gov- 
erned by shear flow accompanied by viscous diffusion, 
which results in the algebraic decay obeying t~^l'^^^. 
However, in a fluid confined by walls with stick boundary 
conditions, both the viscous diffusion and sound propa- 
gation are affected by the walls. The viscous diffusion 
is hampered at the walls, and the accompanying shear 
flow gradually diminishes; however, the sound wave is 
repeatedly reflected at the walls and spreads diffusely^ii^. 
Consequently, the particle motion is finally governed by 
fiow accompanied by the spreading sound wave. Espe- 
cially in a fluid confined between two parallel plane walls, 
the particle velocity relaxation presents a negative t~"^ 



long-time decay differently than that presented in a bulk 
fiuid, which is derived by the mode-coupling theory?^ and 
Green's function method'*. 

In the present study, we investigate the dynamics of a 
single spherical particle in a fluid conflned by two parallel 
plane walls with stick boundary conditions using a direct 
numerical simulation approach. We use the smoothed 
profile method (SPM)iiii2,, which is applicable to a com- 
pressible fluidi^ii^. In SPM, rigid fixed wall boundaries 
can be imposed similar to the representation of rigid par- 
ticleslS. The accuracy of SPM for the present system is 
confirmed by calculating the steady-state mobility of the 
particle, which is compared with the approximate analyt- 
ical solutions. We examine the velocity relaxation of the 
particle after adding an impulsive force in the direction 
parallel to the walls. The velocity relaxation function 
corresponds to the velocity autocorrelation function in 
a fluctuating system. We first consider an incompress- 
ible fluid to investigate the wall effects on the dynamics 
governed only by viscous diffusion. We then consider 
a compressible fluid and investigate the particle motion 
affected by the spreading sound wave arising from reflec- 
tions at the walls. 



II. MODEL 

We consider a system in which a single particle is dis- 
persed in a Newtonian fluid conflned by two parallel plane 
walls, as described in Fig. [1] The stick boundary condi- 
tions are imposed on the walls. The motion of the particle 
is governed by Newton's and Euler's equations of motion 
as 
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M—V = 
dt 



dt 



dl^=^' 



n = N" + N 
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where R, V, and fl are the position, translational veloc- 
ity, and rotational velocity of the particle, respectively. 
The particle has a mass M and a moment of inertia /. 
A hydrodynamic force and a torque are exerted 
on the particle by the fluid, and a repulsive force 
prevents the particle from overlapping with the walls. A 
force F^ and a torque 7V^ are externally applied. The 
hydrodynamic force and torque are evaluated by simul- 
taneously considering the fluid motion. 

The dynamics of fluid are governed by the following 
hydrodynamic equations: 



dp 

dt 



^ + V • {pvv) 



V • {pv) = 0, 



opfp 



(3) 
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where p{r, t) and v(r, t) are the mass density and velocity 
fields of the fluid, respectively. The stress tensor is given 
by 



-pi + r][Vv + [Vvf] 



Vv 



g?? ) (V • V)I, 



(5) 
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FIG. 1. Geometry of the present model system. A single 
spherical particle is in a fluid confined by two parallel plane 
waUs. The particle radius is a, and the wall spacing is H. 
The particle is located at a height I from the lower wall. The 
system is considered to be in a simulation box of size Lx x 
Ly X Lz. The periodic boundary conditions are imposed on 
the end of the box in all directions. The walls are set on the 
top and bottom of the box in the z direction with thickness 
{Lz — H)/2. The walls generate anisotropy for two directions: 
parallel and perpendicular to the walls. In this figure, the x 
and y directions are the degenerate parallel directions, and 
the z direction is the perpendicular direction. 



where p(r, t) is the pressure, 77 is the shear viscosity, and 
r]y is the bulk viscosity. The body force p(j)pfp is added 
such that the rigidity of the particles is satisfied. The 
external force field pcj^wtw is also introduced to impose 
the stick boundary conditions by two parallel plane walls 
such that the force pipwfw constrains the velocity field 
in the wall region to be zero. Additionally, we assume a 
barotropic fluid described hy p = p{p), with the constant 
speed of sound c being 



dp 



(6) 



Equations ©-(H]) are closed for the variables p, v, and p 
without consideration of energy conservation. 

We use the SPM for the direct numerical simulations 
in the present study. The system is composed of three 
regions: the fluid, particle, and wall. In SPM, the 
boundaries between the fluid region and the other re- 
gions are expressed by the continuous phase proflle func- 
tion. There are no boundaries between the particle and 
wall regions because the particle does not penetrate into 
the wall. For this purpose, we introduce a smoothed pro- 
file function (f>x{r,t) e [0, 1], where the index X signifies 
the region of the particle P or wall W. The function (f>x 
represents the boundary between the regions of the fluid 
and A, such that (j)x = 1 in region X and (jjx = in 
the other regions. Using the smoothed profile function, 
the regions of the fluid and X are smoothly connected 
through a thin interfacial region with thickness ^. The 
detailed mathematical expression of 4>x is given in a pre- 
vious paper—. 



III. NUMERICAL RESULTS 

Numerical simulations are performed for a three- 
dimensional box with periodic boundary conditions. The 
space is divided into meshes of length A, which is a unit 
length. The units of the other physical quantities are 
defined by combining 77 = 1 and po = 1 with A, where 
Po is the fluid mass density at equilibrium. The other 
parameters are set to a = 4, ^ = 2, pj, = 1, ?/„ — 0, and 
h = 0.05, where pp is the particle mass density, and h 
is the time increment of a single simulation step. The 
geometry of the present system is described in Fig. [1] 
The particle is set on the midway between two walls as 
I = H/2. The ratio of the wall spacing to the particle 
radius, H* = H/a, is set to various values. The system 
size is Lx X Ly X Lz = 512 x 512 x 32 for H* < 8 and 
Lz = 64 for H* > 8. 

We investigate the relaxation of the particle velocity 
after exerting an impulsive force at the center of the 
particle. The impulsive force is assumed to be suffi- 
ciently small such that the Reynolds and Mach numbers 
of the fiow are sufficiently low. We set the impulsive 
force to produce an initial particle Reynolds number of 
Rcp = 10"'^. In the considered system, two directions 
are specified: parallel and perpendicular to the walls; 
therefore, assuming a low Reynolds number fiow, the re- 
laxation of the particle velocity is described as 



(7) 



(8) 
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where P is the impulsive force exerted on the particle at 
t = and z is the unit vector in the z direction. The 
velocity relaxation tensor ^(i) also depends on the wall 
spacing H and the distance of the particle from the lower 
wall I. In the present study, we focus on the parallel mo- 
tion of the particle 7"(i); therefore, in the following sec- 
tion, we represent the velocity relaxation function in the 
parallel direction by "f{t) for simplicity. For this reason, 
collisions of the particle against the wall will not occur, 
and the direct interactions between a particle and wall, 
including the overlap repulsion force, are not considered 
in the present simulations. According to the fluctuation- 
dissipation theorem, the velocity relaxation function is 
equivalent to the velocity autocorrelation function in a 
fluctuation system: 



lit) 



M 



(9) 



where ks is the Boltzmann constant, and T is the ther- 
modynamic temperature. 

The important time scales in the dynamics of a single 
particle are those of viscous diffusion and sound propa- 
gation over the length of the particle radius: Ti, = a} jv 
and Tc — a/c, respectively, where v = tj/pq is the kine- 
matic viscosity. We define the compressibility factor by 
the ratio of these time scales as 



T,, 



1/ 

ac 



(10) 



which provides a measure of the importance of sound 
propagation in the dynamics of a single particle. In the 
present simulations, we adjust the fluid compressibility 
by the compressibility factor. 



A. Steady-state mobility in an incompressible fluid 

First, we estimate the particle mobility /i in steady 
state for accuracy testing. In these simulations, a con- 
stant force was continuously exerted on the particle in 
the direction parallel to the walls. The mobility was cal- 
culated as the ratio of the particle velocity to the added 
force after reaching steady state. The mobilities in such 
situations have been studied by Faxen using analytical 
theoriesiii^. The solution was given as the power series 
of the ratio of the wall spacing to particle radius, in which 
the first few terms were derived as 

6mjan{H) = 1 - 1.004A + 0.418A^ 

+0.21X^~0.169X^ + 0{\'^), (11) 

where A — 2/H*. The simulation results are presented 
in Fig. [5] and are compared with the approximate solu- 
tions given by Eq. (|11|) . The simulations were performed 
with the various system side lengths in parallel direc- 
tions to the walls, and Ly [L^ — Ly), to investigate 
the system size effects due to the periodic boundary con- 
dition, which are diminished to an insignificant level at 
Lx > 256. However, the simulation results underestimate 
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FIG. 2. Wall spacing dependence of the mobility of a single 
particle in a fluid conflned by two parallel plane walls. The 
length of the simulation box sides parallel to the walls are 
changed to ~ Ly — 64, 128, 256, and 512 to investigate 
the effect of the periodic boundary condition. The broken line 
represents the analytical solution given by Eq. (|lip . 



the mobility even at — 512; the deviations from the 
solutions given by Eq. ([TT|) are less than 2% for iJ* > 6 
and approximately 5% only for H* = 4. As for H* ~ 4, 
however, because the order of the power A^ is still 1, 
a truncation error of a few percent arises in Eq. (|11[) . 
Therefore, the simulation error for H* = 4 can be less 
than 5%. 

From the present results, the walls with stick bound- 
ary conditions are successfully implemented by SPM. The 
implementation of confinement by different geometries is 
also available in SPM. 



B. Velocity relaxation in an incompressible fluid 

Here, we consider an incompressible fluid as a solvent 
fluid. To treat an incompressible fluid, we assume an 
infinite speed of sound in Eq. ([6|) and ignore the deviation 
of the fluid density. With this assumption and the mass 
conservation law Eq. ([3]), the solenoidal condition for the 
velocity field is derived as 



V ■v = Q. 



(12) 



Therefore, Eqs. (HJ, ([5]), and (fT2|) were solved in combi- 
nation as the hydrodynamic equations. 

The simulation results of the velocity relaxation func- 
tions for wall spacings H* = A and 6 are presented in 
Fig. [21 The relaxation functions decrease monotonically, 
and the reduction rate increases with a decrease in the 
wall spacing H*. Exponential long-time decay is ob- 
served, as derived from the analytical theories^^. Such 
decay contrasts the relaxation function in a bulk fiuid. 
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FIG. 3. Velocity relaxation functions of a particle in an in- 
compressible fluid conflned by two parallel plane walls. The 
wall spacings are H* — 4 and 6. The analytical solution of 
the corresponding function in a bulk fluid is also presented^^. 
The bold dotted line represents the long-time tail At~^^^, as 
given by Eq. (fT5)) . 



which describes long-time decay obeying the power law 
given bj^ 

^buik(^)^ 1 Pp/M^/^ as t^oo. (13) 
9y/n po \ t / 

The exponential decay of the relaxation function is ex- 
pected to reflect the loss of fluid momentum at the walls 
with stick boundary conditions. Therefore, the time 
when the relaxation function in a confined fluid starts to 
remarkably fall below that in a bulk fluid is related to the 
time scale when the fluid flow generated by the particle 
reaches the walls. In an incompressible fluid, the tem- 
poral evolution of the fluid flow is only accompanied by 
viscous diffusion, whose time scale over the distance be- 
tween the particle surface and the wall {H/2 — a) is given 
by tI = {H/2 - afjv = (H* - ifr^. For the wah 
spacings H* — 4 and 6, the viscous diffusion time scales 
are tI/t^ = 1 and 4, respectively. As demonstrated in 
Fig. [31 the starting time of the deviation of the function 
from that in a bulk fluid corresponds to the time scales 
' 1/ ■ 

The velocity fields around the particle with various 
wall spacings at the time t/r^ — 2.81 are displayed in 
Fig. 01 Compared with a bulk fluid, the attenuation of 
fluid velocity is clearly observed in a confined fluid with 
H* — A, corresponding to the decay of the shear flow or 
the vorticity by the walls. However, the attenuation of 
fluid velocity is not observed for H* — 6. This result co- 
incides with the fact that the time t/r^ = 2.81 is earlier 
than the viscous diffusion time scale for H* = 6, namely, 
rl/Ty = 4. The walls restrict the viscous diffusion of the 
flow field to form a laminar flow, which corresponds to 
the extended distribution of the vorticity along the walls. 
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(b)H =6 
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FIG. 4. Velocity fields around the particle at t/T^ = 2.81 
in confined fluids with a wall spacing of (a) H* = 4 and (b) 
H* = 6 and (c) in a bulk fluid are shown. The walls are set on 
the upper and lower boundaries of the images for the confined 
fluids, and the corresponding region is presented for the bulk 
fluid. The cross-sections presented here are perpendicular to 
the walls and parallel to the impulsive force direction and in- 
clude the particle center. The direction of the impulsive force 
is to the right in the images, and the particle is represented 
by a black circle. The vorticity of the velocity field V x d is 
described using a color scale, moving from negative (darker) 
to positive (lighter) vorticity. The vorticity is normalized by 
the factor of r^/Rcp. 

C. Velocity relaxation in a compressible fluid 

In a compressible fluid, the velocity relaxation func- 
tion exhibits essentially different behaviors from those in 
an incompressible fluid. The simulation results of the re- 
laxation functions are presented in Fig. [5] Remarkable 
differences are the sign inversion and the subsequent neg- 
ative long-time algebraic decay. The power of the long- 
time decay is and the coefficient is analytically de- 
rived as^i^ 

In the present simulation results, the asymptotic ap- 
proaches to Eq. (HH) are accurately reproduced. In a 
bulk fluid, the long-time decay is positive and propor- 
tional to t~^/'^, given by Eq. (|13p as an incompressible 
fluid. Although the sign inversion can occur for a large 
compressibility factor as e — 1.5, the relaxation function 
eventually becomes positive. 

For the compressibility factor e — 0.6, the temporal 
evolution of the velocity field in a confined fiuid is com- 
pared with that in a bulk fluid in Fig. [51 The flow ac- 
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FIG. 5. Velocity relaxation functions of a particle in a com- 
pressible fluid confined by two parallel plane walls with the 
wall spacing H* = 4. The compressibility factors are (a) 
e — 0.6 and (b) e = 1.5. The analytical solutions of the cor- 
responding functions in a bulk fluid are also presente d^^''^° . 
The functions are described by solid lines (positive values) 
and dashed-dotted lines (negative values). The bold dotted 
lines represent the positive long-time tail At~^^^, given by 
Eq. (|13|) . The bold dashed-two dotted lines represent the 
negative long-time tail —Bt~^, given by Eq. (|14[) . 



companied by a sound wave is source-sink flow: the flow 
moves from the source (positive divergence region) to the 
sink (negative divergence region). The divergence of the 
velocity field corresponds to the fluid density deviation 
as 



Dt po 



(15) 



In a bulk fluid, a sound wave propagates away from the 
particle; therefore, eventually, the effect of the sound 
wave on the particle motion disappears, and the particle 
motion is governed by shear flow accompanied by viscous 
diffusion. However, in a confined fluid, a sound wave is 
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FIG. 6. Temporal evolution of the velocity fields in a con- 
fined fluid with H* = 4 (upside in each time) and in a bulk 
fluid (downside in each time). The compressibility factor is 
e = 0.6. The walls are set on the upper and lower boundaries 
of the images for the confined fluids, and the corresponding 
regions are presented for the bulk fluids. The cross-sections 
presented here are perpendicular to the walls and parallel to 
the impulsive force direction and include the particle center. 
The direction of the impulsive force is to the right in the im- 
ages, and the particle is represented by a black circle. The 
divergence of the velocity field V • i) is described using a color 
scale, moving from negative (darker) to positive (lighter) vor- 
ticity. The divergence is normalized by the factor of r^/Kep. 



repeatedly reflected at the walls and continuously affects 
the particle motion. Multiple reflections at the walls re- 
sults in the spreading of sound wave, which is known as 
the overdamped diffusive modeS-^i. Because viscous dif- 
fusion is prevented by the walls and the corresponding 
shear flow gradually disappears, the particle motion is 
eventually governed by the sound wave. The spreading 
sound wave is associated with backward fluid flow and 
causes backtracking of the particle, which corresponds to 
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FIG. 7. Velocity relaxation function of a particle in a compressible fluid confined between two flat walls. The wall spacings are 
(a, b) H* = 4 and (c, d) H* = 6. The compressibility factor has values of e = 0.3, 0.6, 1.0, and 1.5. The simulation results 
are represented by solid lines (positive values) and dashed-dotted lines (negative values). The bold dashed- two dotted lines 
represent the negative long-time tail —Bt~^, as given by Eq. (|14[) . 



the negative velocity relaxation function. 

The long-time decay of the particle velocity reflects the 
eventual mechanism transmitting the fluid flow: viscous 
diffusion or sound propagation. In a bulk fluid, the parti- 
cle velocity finally decays, obeying Eq. (|T3)) . which corre- 
sponds to viscous diffusion of the shear flow. The volume 
of the viscous diffusion region at time t is proportional 
to (j/<)'^/^. In a confined fluid, the different long-time 
decay Eq. ([T4| appears. This negative decay originates 
from the spreading sound wave, whose spreading volume 
at time t is proportional to H{ct)^. 

In Fig. [3 the relaxation functions with wall spacings 
H* ^ A and 6 for various compressibility factors are dis- 
played. Discrepancies of long-time decay with Eq. 
for small compressibility factors are presumed to result 
from system size effects due to the periodic boundary 
conditions, which will appear after the time when the 
sound wave generated by the particle reaches the end 



of the system: t = {Lx/2 — a)/c — {Lx/2a — l)eT^. 
The time when the sign inversion occurs changes non- 
monotonically with respect to the fluid compressibility. 
With an increase in the compressibility factor, the time 
of the sign inversion occurs earlier when the compress- 
ibility factor is small; however, the sign inversion occurs 
later when the compressibility factor is large. In other 
words, there is a compressibility factor at which the time 
of the sign inversion is minimized, which we call the char- 
acteristic compressibility factor. The characteristic com- 
pressibility factor depends on the wall spacing H*; for 
example, from Fig. [71 the characteristic compressibility 
factors are presumed to be e ~ 0.6 for H* = 4 and e sa 1.0 
for H* = 6. More detailed evaluations of the character- 
istic compressibility factors for various wall spacings is 
provided in Fig. HI 

The sign inversion of the velocity relaxation function, 
namely, the backtracking of the particle, is expected to 
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FIG. 8. Classification of calculated velocity relaxation func- 
tions for various wall spacings and compressibility factors. 
The time when the sign inversion of the velocity relajcation 
function occurs non-monotonically changes with respect to 
the compressibility factor. For each wall spacing, with the in- 
crease in the compressibility factor, the decrease and increase 
of the time of the sign inversion are represented by upward and 
downward triangles, respectively. The time of the sign inver- 
sion is minimized at the characteristic compressibility factor, 
which is represented by open diamond. The broken line rep- 
resents Eq. (|16|) with = 1. The dashed double-dotted line 
represents e = 1.34. 



occur after the sound wave reflected at the wall reaches 
the particle. The time scale of this event is estimated by 
t] = 2{H/2 - a)/c= {H* - 2)eT^. Moreover, sufficient 
decay of the shear flow around the particle is required, 
and the time scale of the shear flow decay is t}, , as given 
in the previous section. Here, we define the confined 
compressibility factor as the ratio of the two time scales: 



4e 



H* -2 



(16) 



When the confined compressibility factor is small, e'l' < 1, 
the spreading of the sound wave caused by reflection at 
the walls progresses faster than the decay of the shear 
flow. For smaller confined compressibility factors, the 
sound wave spreads and is weakened more rapidly; there- 
fore, further reduction of the shear flow is required to 
cause backtracking of the particle. In short, when e^^ < 1 
is satisfled, backtracking of the particle occurs at a later 
time for a smaller confined compressibility factor. How- 
ever, when the confined compressibility factor is large, 
> 1, the reflected sound wave reaches the particle af- 
ter sufficient decay of the shear flow; therefore, particle 
backtracking occurs at a later time for a larger confined 
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FIG. 9. Velocity fields around the particle at t/r^ = 2.81 (a) 
in a confined fiuid with wall spacing H* = 4 and (b) in a bulk 
fiuid. The compressibility factor is e = 1.5. The walls are set 
on the upper and lower boundaries of picture for the confined 
fiuid, and the corresponding region is exhibited for the bulk 
fiuid. The cross-sections shown here are perpendicular to the 
walls and parallel to the impulsive force direction with includ- 
ing the particle center. The direction of the impulsive force 
is to the right in the pictures and the particle is represented 
by a black circle. The divergence of the velocity field V • v 
is described by a color scale, going from negative (darker) to 
positive (lighter) vorticity. The divergence is normalized by 
the factor of r,//Rep. 



compressibility factor. Consequently, the characteristic 
compressibility factor is expected to satisfy e'^ ^ 1. In 
Fig. [51 the line of Eq. with = 1 almost coincides 
with the characteristic compressibility factors at H* < 6; 
however, this relation fails at H* > 6, where another 
mechanism of backtracking must be considered. 

For the compressibility factor e = 1.5, a comparison 
of the velocity fields in confined and bulk fiuids at time 
t/r^ — 2.81 is presented in Fig. [3] Although this time 
is earlier than the time rj, the relaxation function in 
the confined fiuid is negative, as illustrated in Fig. [5l^b); 
therefore, the backtracking is caused by the pressure from 
the sound wave remaining in the vicinity of the particle. 
Backtracking by this mechanism is also observed in a 
bulk fluid when the compressibility factor is sufficiently 
large22i^, and the condition for which the backtracking 
occurs is estimated as e > 1.34 from the analytical so- 
lution— i^. There is only a slight difference in the di- 
vergence of the velocity field between the confined and 
bulk fiuids; however, in a confined fiuid, because shear 
fiow is diminished by the walls, backtracking by a sound 
wave can occur at an earlier time than in a bulk fiuid, as 
demonstrated in Fig. [5l^b). Considering this mechanism 
of particle backtracking, the characteristic compressibil- 
ity factor will eventually converges to e = 1.34 with an 
increase in the wall spacing H* . The successive change in 
the dependence of the characteristic compressibility fac- 
tor on the wall spacing, namely, from = 1 to £ = 1.34, 
is observed in Fig. [S] The wall spacing at which this 
crossover occurs is estimated as H* sa 7.4. 
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IV. CONCLUSION 

We investigated the dynamics of a single particle in 
a fluid confined by two parallel plane walls using SPM 
to perform direct numerical simulations. In particular, 
we calculated the velocity relaxation of the particle after 
an impulsive force was added in the direction parallel to 
the walls. The velocity relaxation function corresponds 
to the velocity autocorrelation function in a fluctuating 
system. 

In an incompressible fluid, the relaxation function de- 
cayed more rapidly in a confined fluid than in a bulk fluid; 
the long-time decay was exponential differently from the 
power law t""^/^ observed in a bulk fluid. The rapid de- 
cay reflects the loss of shear flow due to the hindrance of 
viscous diffusion by the walls. Therefore, the time when 
the relaxation function in a confined fiuid falls below that 
in a bulk fluid corresponds to the time scale of viscous 
diffusion over the distance between the particle surface 
and the walls, t^. 

In a compressible fluid, sign inversion and subsequent 
negative long-time decay obeying were observed 
in the velocity relaxation functions. A corresponding 
spreading sound wave arising from the multiple reflec- 
tions at the walls was also observed. The time of the sign 
inversion changed non-monotonically with respect to the 
compressibility factor; that is, the time of the sign inver- 
sion was minimized at the characteristic compressibility 
factor. The characteristic compressibility factor satisfles 
e'' — 1 when the wall spacing is small and is given by 
e = 1.34 when the wall spacing is large, where the con- 
fined compressibility factor is defined by Eq. ((T5)) . The 
crossover of the dependence of the characteristic com- 
pressibility factor on the wall spacing is estimated to be 
H* « 7.4. The backtracking of the particle is caused by 
the spreading sound wave arising from reflections at the 



walls when the characteristic compressibility factor satis- 
fies = 1 and is caused by the remaining sound wave in 
the vicinity of the particle when the characteristic com- 
pressibility factor is given by £ = 1.34. Backtracking via 
the latter mechanism can also occur in a high compress- 
ible bulk fluid. 
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